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Abstract 

Nowaday, in study of effective interactions, more attention is devoted to single-particle properties 
of near-magic nuclei and bulk properties of deformed ones but quasiparticle states of the latter 
are rarely used so far because of theoretical difficulties. In particular, the angular momentum 
projection remains too time-consuming for such calculations and the methods, which are based on 
the transformation to an intrinsic frame have some unsolved problems such, e. g., as quantum 
fluctuations of rotational recoil in the description of quasiparticle-rotation coupling. To remove a 
part of these difficulties, the method of the optimal intrinsic frame-of-reference is developed. After 
applying the Mikhajlov transformation to obtain the nuclear Hamiltonian in the intrinsic frame, 
approximate constraints on nucleon's variables are substantiated, and the quasiparticle structure of 
the nucleus orientation angle operators is investigated. That gives possibility to use the variational 
principle to derive equations for matrix elements of these operators. An approximation similar to 
the cranking model (CM), but with the quantum rotational recoil, is formulated, which may be 
considered as a generalization of the usual self-consistent CM. Simplified model calculations for 
rotational bands in ^^^Er show that taking into account the recoil operator considerably improves 
the agreement with experimental data. 

PACS numbers: 21.60.Ev, 21.60.Jz 
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I. INTRODUCTION 



When studying nuclei far away from the stabihty valey, a need arises to improve the meth- 
ods of the microscopic description of deformed nuclei since, on the one hand, many nuclei in 
the region are deformed, and on the other hand, the semi-microscopic methods using phe- 
nomenological single-particle potentials and residual interaction are often inapplicable here 
because of lack of data for determining the model parameters. The microscopic description 
of traditional regions of deformed nuclei, near the stability valey, has a long history (see, 
e. g., books jll - Pi, reviews and references therein) with the latest time witnessing 

great progress by virtue of improvements in computer methods. However, even for these 
regions, microscopic methods are applicable with considerable restrictions; in addition, a 
amount of self-consistent microscopic calculations for heavy nuclei is small so far, especially 
as compared with the huge amount of existing experimental material. The development 
of microscopic methods for such nuclei is of interest for many reasons, and in particular, 
because this opens up new possibilities for studying the nuclear effective interactions. 

At present in study of effective interactions and even in ab-initio calculations, the apparent 
shift of attention from ground-state bulk to single-particle properties of magic and near- 
magic nuclei occurs (see, e. g., [l^-[l8|). At the same time, when determining interaction 

amnd states, just the deformed nuclei rather than magic 
, but the quasiparticle states of the deformed nuclei are 
rarely used so far because of theoretical difficulties. 

In the deformed nuclei, particle- vibration coupling is much weaker, in many cases, than in 
spherical ones as the main part of this coupling is included in the mean deformed field l| (p. 
242). Therefore, along with vibrational states, there are many excited states of essentially 
pure single-quasiparticle nature in deformed nuclei (the probability of phonon admixtures 
are less than 5% j48|; an evaluation of amplitudes is given, also, in Ref. jl]), that makes, 
to some extent, these nuclei similar to the magic ones. The treatment of the superfiuidity, 
which gives single-quasiparticle spectra in deformed nuclei (instead of single-particle ones in 
the magic nuclei), is elaborated enough well. But in addition, there exists the quasiparticle- 
rotation interaction in the deformed nuclei. The matrix elements (ME) of this interaction 
depend on the effective forces much weaker than the ME of the quasiparticle-phonon in- 
teraction in spherical nuclei, therefore the infiuence of the forces is more transparent. The 
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main difficulties in theoretical description of these nuclei are related to the restoration of the 
broken rotational symmetry, that is a serious obstacle in self-consistent microscopic calcula- 
tions. Therefore the overcoming these difficulties would have important advantages. Then 
the microscopic description of deformed nuclei and especially of odd ones could yield vast 
valuable information on effective interactions: degree of their universality, their dependence 
on nuclear shape and excitation. 

For the microscopic description of deformed nuclei, there are used two ways: the ffist 
method exploit wave functions with broken symmetry (as a leading approximation or as 
some building elements), the second one works with wave functions which on every stage 
have symmetries of the nucleus Hamiltonian, are used. The ways are described in many good 
reviews, only some of them are pointed out here: {4]- isl. The brief listing^f methods with a 
description of their advantages and drawbacks may be found, e.g., in Ref. [9|. Supplementing 
these reviews it worth mentioning the self-consistent collective coordinate method and other 



approaches within the framework of time-dependent Hartree-Fock (HF) method (see 2ll |22| 



and references therein) and also the interacting boson model, combined with the cranking 



one, [2J|, pseudo-SU(3) model [2J|. Some of publications, more closely related to the given 
investigation, are mentions below in the given paper. 

Since deformed potentials enable to allow for an essential part of multinucleon correlations 
in a simple form, the approaches that use the symmetry violation prove to be most efficient. 
Among the methods using such approaches, the angular momentum projection (AMP) and 
the intrinsic frame-of-reference (IF) are basic sj. Both AMP and IF methods bear on some 
common physical ideas, in particular, the concept of intrinsic wave functions and deformed 
potentials in IF (the projection of wave functions can be interpreted as a transformation 
of the functions to IF, complimented with integration over all possible orientations of the 
frame). Nevertheless, mathematical formulations of the methods differ markedly. To date, 
AMP is much more developed and applied than IF. The development of computers has made 
possible self-consistent calculations by configuration mixing of states obtained by variation 
after AMP. Only some publications of such works performed in the frame work of both the 



nonrelativistic and relativistic approaches, will be pointed out here: [25|-[30||. Nonetheless, 
further development of microscopic theory of deformed nuclei is needed, since the application 
of self-consistent AMP method is so far possible only with considerable restrictions, espe- 
cially it relates to describing the interaction of single-quasiparticle and rotational degrees of 



freedom in heavy nuclei. 

The second of the mentioned approaches using the symmetry violation, the IF method, 
finds practically no application in self-consistent microscopic calculations. Moreover, despite 
numerous studies (see, e. g., jst, jsih 33| and references therein) and the apparent, at first 



sight, physical obviousness of the method, the description of a nucleus in IF encounters still 
with serious unsolved problems which attenuate the activity in the direction. As a result, 
the basic methods for heavy odd nuclei are still semi-classic or semi-microscopic ones. These 
are the cranking model (CM) and its more general version, tilted- axis CM (see reviews 0,0] 



and latest papers 



3j]-|37|) and also a many-particle plus rotor model (e. g., |38l. |39|). 

In the present paper, an attempt is made to overcome some of the difficulties of the 
IF method. This can be useful for the formulation of an efficient approach to describe 
the quasiparticle-rotation interaction in nuclei (one of advantages of IF method is that it, 
unlike the AMP, requires no calculation of overlap integrals), for better understanding of the 
relationship between the AMP and IF, and for the derivating quantum corrections to widely 
used models, first of all, to the CM. The neccessity for such corrections calls for a special 
explanation. 

In the low approximation the AMP and IF methods may be reduced to the self-consistent 
CM (SCM). Within the framework of the first method, the SCM equations are derived on the 



basis of the Kamlah expansion 40| for the expectation value of the nuclear Hamiltonian over 



the projected wave functions. A series arising here descends in the degree of coherence of its 
terms [3| . The second method gives as a series not for the expectation value but the nuclear 



Hamiltonian, transformed to variables of rotating frame of reference [33|. Both methods 
derive the SCM equations with taking account of only two first terms of the corresponding 
series. One of the terms describes, conventionally speaking, the intrinsic motion, since it is 
independent of rotational variables, and the second one describes the quasiparticle-rotation 
interaction. Conventionality of such terminology consists in the fact that the expectation 
value of the Hamiltonian is found by wave functions, allowing for rotation; moreover, the 
effective nucleonic forces may depend on nuclear density, which itself may depend on an- 
gular momentum of the nucleus. The rest of the terms, ignored in the SCM, describe the 



rotational recoil effect, more precisely, its fluctuation part, since some averaged part of t 



he 



effect is accounted for in the quasiparticle-rotation interaction of the SCM (see, e. g., 4l[|). 



In the other terminology, the rotational recoil effect, like a similar effect for translations. 



corresponds to subtraction of the spurious motion in IF. 

Importance of the quantum recoil effect has been discussed many times (see, e. g., 



Refs. 



41 



43|): yet, majority of works that use the SCM ignores the effect. Corrections to 



the SCM, derived on the AMP basis in Ref. |44|, are applicable only to states with large val- 
ues of the quantum number K, i. e., the projection of the angular momentum of the nucleus 
to the symmetry axis. An efficient means to correct the SCM on the AMP basis is the pro- 
jecting of cranked states (see Refs. [36|, l37|, |45| and references therein). Here, the projection 



before variation, necessary to reproduce the correct moment of inertia, is replaced by its first 
approximation, SCM, and the projection (after variation) of obtained functions gives quan- 
tum corrections to the SCM. However, such approach can so far be implemented only in the 
HF approximation without pairing correlations and does not take into account an influence 
of the quantum rotational recoil on mean fields, which is essential in both particle-hole and 
particle-particle channels. 

For odd nuclei, probably, for the present there is no publication of a calculation, allowing 
for the quantum effect of rotational recoil on the basis of the SCM or a generalization of it. 



At the same time, as calculations have shown 



46| . the SCM yields some underestimation of 



quasiparticle-rotation coupling as compared with the observed one. The difference between 



the model used in 



46j and the SCM is immaterial for this conclusion (the effect of the 



perturbative treating of 3-quasiparticles states in the framework of SCM was considered in 
Ref. (47|). The underestimation is quite opposite, though less by magnitude, to the well- 
known overestimation of the Coriolis coupling in the particle- rotor model It shows up in 
theoretical rotational level energies for strongly-mixed bands which, in the SCM, exceed the 
experimental ones with the differences increasing for the increasing nucleus spin /. Note that, 
at first sight, the application of the SCM in the independent-quasiparticle approximation to 
description of positive parity band in I55,i59,i6i£)y leads to agreement with experiment, which 
has been demonstrated in Ref. |47j], and was confirmed by our calculations. However, the 
agreement arises due to ignoring another, important for these states, effect. The matter is 
that, in I55,i59,i6i£)y^ ^j^g states of intrinsic motion, originating from shell ^13/2, contain large 
quasiparticle- vibrational admixtures ( > 20% in states 3/2+ [651] and 7/2"*" [633] |48| ). When 
the admixtures are accounted for, the ME of the quasiparticle-rotation interaction decrease 
(e. g., for ^^^Dy, such attenuation of the interaction was discussed in the framework of the 
particle- rotor model |4^) and the agreement of the SCM with experiment is broken. 



Let us return to the description of goals and features of the present work. The work 
develops the method of optimal IF (OIF), based on the Mikhajlov [? ] transformation of 
wave functions and operators from the laboratory coordinate system to a rotating one (see 



5Ch52|). To this effect, approximate conditions of constraint on nucleonic variables are sub- 



stantiated, the quasiparticle structure of the nucleus orientation angle operators is studied, 
the equations for ME of these operators are derived, and an approximation is stated, which 
can be considered as a CM generalization, containing the quantum effect of the rotational 
recoil. The results obtained are intended for odd nuclei, though most part of them can be 
useful for even ones, too. 

The work sets no goal to carry out microscopic calculations for specific nuclei, as this is 
a separate serious problem. The calculations presented in Section [IV] pursue another aim of 
no less importance: to analyze main features of the suggested approach and reveal the role 
of different effects, starting with a simple Hamiltonian and semi-microscopic approximation, 
which combines diagonalization with perturbation theory. The example ^^■^Er is used to show 
that the accounting for rotational recoil operator within the framework of the generalized 
CM improves noticeably the agreement with experiment. 

The paper is organized as follows. In Sect. [Tll after a short description of the transfor- 
mation to the IF, the constrains between nucleonic variables in this frame are considered. 
Then, the general ideas of the OIF method are introduced. The quasiparticle structure of 
the nucleus orientation angle operators is studied in Sect. Illli At the end of this section, the 
difference between the OIF and the approximate AMP methods is shown. The formulation 
of the method in the simple model and the results of calculations in the model are given 
in Sect. HVl Section |V] contains a summary and an outlook to future development of the 
method. 



II. TRANSFORMATION TO THE INTRINSIC SYSTEM AND CONSTRAINS 
BETWEEN THE VARIABLES 



General ideas of the method of a rotating frame of reference are set forth in Ref. [33 1 
where also references to previous works by other authors are given. The method employs 
V. Mikhajlov concept of the unitary transformation of wave functions and operators from 



the laboratory frame to some intrinsic one that rotates with the nucleus (see 50l-l52[| and 



references therein). The principle of the method is in brief as follows. Let Xj, x'^ and 
'^i{xi), denote, respectively, nucleon variables and wave functions of the nucleus 

in the laboratory and intrinsic systems, where i = 1, ... A, A is the number of nucleons, 
/ is the nuclear spin. We denote orientation angles of the intrinsic with respect to the 
laboratory frame as "^a with understanding to be the entire set of and Xi to be the set of 
variables {xi}. The transformation U is defined with the equation: \E'/(xj,C) = "^lix'i,^) = 
U'^ i{x^, d). The additional variables C,a are introduced for the number of variables before the 
transformation to be equal to that after it and can be interpreted as the orientation angles 
of the rotating frame of reference with respect to the axes of the nucleus. The redundant 
degrees of freedom are eliminated with constraints on variables in the IF. 

The orientation of the nucleus with respect to the rotating frame of reference is described 
in terms of x[ by three commuting operators Oa{x'i, ...x'^) = 9a{x[), which are dependent on 
momentum, spin and isospin operators of nucleons as well as on coordinates. The selection of 
operators 9 is carried out from considerations of the problem simplification and is restricted 
by the requirement that their commutation rules with the projections of the nuclear angular 
momentum Ja{x'i) be the same as those for the orientation angles of a rigid rotor da with 
operators la of the projections of its angular momentum to the laboratory axes: 

[/a, = babid), [Ja, ^0b] = bab{0), 9,] = . (l) 

If are Euler's angles or three angular parameters that fix the rotation axis and the rotation 



angle of a rigid body, the expressions for hab{'d) can be found in Ref. 53|. Below, the 



Cartesian components of the vector, defined by the rotation axis and angle, are used as 



'da {a,b = x,y,z); the corresponding expressions for hab{'&) are given in Ref. |33|. The 
transformation U is chosen so that the equations hold: 

J„(xi) = UJ^{x[) U-^ = Ud) = -^K,^. D{C) = UD{d) = D{-9{x[)), (2) 

were D{d) = are Wigner D-functions; the summation over repeating Greek indexes is 

implied here and below. This requirement leads to the expression U = exp {^aJa) exp {—9 pip) 
50|. Note that the use of the operator angles 9a{x'j) presupposes that the nuclear wave func- 
tion before W-transformation, i. e. the function in terms of the laboratory coordinates, has 
the subsidiary variables on the left from nucleonic ones. Then, in the simplest case, the 
function before U/ie[/-transformation has the form \E'/(xj,C~^) = f{C~^)^ and after it 



— = f{9) ^j{x'^,i)). To be more precise, this is valid in the case of small (, which 

is discussed below, while a more involved relation, based on Eq. ([T]), should be used for 
arbitrary (. 



In the 3D case, the nuclear Hami 
by Mikhajlov as an infinite series 



51 



tonian in the rotating frame of reference was obtained 



52|: 



H = UHU 



-1 



y- 



[H, i9aj],i0a2], ■ ■ •] 



(3) 



where = 1^ is the projection of the angular momentum of a rigid body to the 

intrinsic axes, Sa = b'^f^i^) (I/s — J/s), and b^J^(9) is the matrix inverse with respect to bai3{9). 
The series as many other expressions in the approach in question can be radically simplified, 
if the ME of operators 9 are small enough. Then the transformed Hamiltonian (jS]) and the 
commutation relations ([1]) for 9 take the form: 



a a 

[Ja,i9b] = Sab a,b = x,y . 



(4) 
(5) 



The present work is confined to consideration of axially symmetrical nuclei, therefore, op- 
erators 9z are not used in the formulas (some deviations from the axial symmetry can be 
allowed for by mixing states with different projections of the angular momentum to the 
sjTiimetry axis). The smallness of 9 is ensured by imposing constraints on variables. If 
the constraints are chosen so that (a were small enough before the W-transformation, the 
equations for D-functions in (|2]) can be replaced with a simpler equation (a = —9a{x'i), 
whereupon 9a{x'j^) will be small too, and the rotating frame of reference will coincide with 
the intrinsic one. The definition of the intrinsic system will be addressed again below, now 
the quantum constraint conditions come into consideration. As the work in Ref. 33|, the 



present one makes use of approximate constraints, but because of their importance, provides 
more rigorous substantiation and shows what results follow from them. 



Correct, but with a need to overcome great difficulties, methods for ta. 
quantum constraints for the nuclear rotation were proposed in Ref. 
therein) on the basis of Faddeev-Popov functional integral and in Ref. 
BRST-symmetry. We will use the method of Ref. 
statement, indicating only suggested changes to it. 



dug into account 
(and references 
on the basis of 



3l| and confine ourselves to its brief 
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In classical mechanics, constraints can be introduced as follows. If the Lagrangian in the 
laboratory system is assumed independent of additional variables (a, there arise constraints 
for momenta, canonically conjugate to these variables: Ua = dL/d(a = 0. The constraints 
are introduced into the Hamiltonian as an additional term with Lagrange factors, while 
equations Ca = can be chosen as gauge conditions. 

In quantum mechanics, when the vacuum transition amplitude for a nucleus is written 
in terms of a path integral, the additional variables in the laboratory frame are introduced 
into the path integral with the multiplier 

1= f l[6{n,)6i(,)DnaDC,, (6) 

ab 

where 6{. . .) is delta function. Computation of the resulting functional integral can be 
replaced (see [3l| and references therein) with solving the eigenvalue problem for the effective 
Hamiltonian 



Ha = lim 

en p 



a a J 



Here H is the nucleus Hamiltonian in the laboratory system, A is an arbitrary constant, 
having no effect on the final expressions; Ua^C) = —id/dCa are conjugate to the angular 
variables. 

In the Hamiltonian -ffgff, the spurious motion, described by the additional variables C^, 
is separated from the real one and represents the harmonic oscillator motion. We will be 
interested only in those eigenstates of the total Hamiltonian ([7]), which correspond to the 
ground state of the ghost motion. At D — )■ 0, the states are separated by infinite energy from 
all states, where the spurious oscillator is excited. After transformation to the IF variables, 
the nucleonic and supplementary variables in the Hamiltonian (and in wave functions) are 
not separated any more, but the limit P — )■ ensures the fixing of the redundant degrees of 
freedom and the absence of spurious motion. For the ground state of the oscillator, we have 



(Ca) = l/'2,^J A/V and, with A fixed, the relationship {Q) oo at — )■ is valid, which 
makes it impossible to ensure smallness of da- Yet, the freedom in selecting A allows to take 
it such that ^ — )■ at — and the ratio A/T) were small enough. This ensures sufficient 
smallness of {Ca)i ^^"^ therefore, after the W-transformation, smallness of {0"^)- Note that 
one condition ^ — suffices for the approximation used below, whereas the limit P — )■ is 

n 

introduced here to trace a parallel with Ref. |31| . 



Since at C — ^ 0, the matrix {(), inverse with respect to bab{C), reduces to Kronecker 6ab 
symbol and since we will use only those states of the Hamiltonian H^q, where the spurious 
oscillator is in the ground state, one can replace i7a(C) -^a(C) iii Eq. ([7]) by the operator 
77/3 . Validity of this change is especially clearly seen in the path integral, 
corresponding to the Hamiltonian ([7]), because of the presence of 5(Ca) in the integral. Con- 
sidering the above and the commutation relation [il/? ,b~^] — )■ at (a — ^ 0, one can 
finally replace the operator IJalla in Eq. ([7]) by the operator lala- Here Ea/^-y is the unit 
antisymmetric tensor and Ia{i^) = is the projection of the angular momentum 
of a rigid body to the intrinsic axes, = hp^II^. 

After the transformation to IF variables {xj, Ca} — > {x'i, '&a\, operators, involved in Eq. ([7]), 
have the form /a(C) = lai.'^) — Jai^'^, . . . ,x'j^), (a = —^a{x[, . . . ,x'j^), and the effective 
Hamiltonian itself is described by the expression 



K a a . 



Note that the effective Hamiltonian in the approach 31| resembles Eq. (jH]), but contains the 
Hamiltonian H instead of H. Therefore the description of a nucleus in IF is for the most part 
accomplished only due to the constraints and the precision of the method is strongly limited 
with accuracy of accounting for constraints. In the suggested approach, as will be seen below. 



good results are produced. 
The method of Ref. 



Dy accounting for constraints even in the lowest approximation. 
31| . in its nowaday formulation, is practically inapplicable to self- 
consistent description of nuclei, as it requires an exact compensation of divergent (in the 
limit P — 7- 0) terms in every order of the perturbation theory. Used below is another, 
more primitive method of accounting for constraints, which nonetheless contains the most 
important effects, described by the Hamiltonian (|H]). Instead of the Hamiltonian H^q, we will 
use H, but with account of two important properties of eigenfunctions \I'/ of the Hamiltonian 



ifgg. Firstly, the expectation value {"^j 
9a should be small, and, secondly, (see 



9' 



46 



y f) should be small, and, hence, ME of operators 



52|) 



(*,|IJ|^,) = /(/ + !) . (9) 



The second of the conditions is simpler, so we start discussing it first. Correctness of ([9]) 
can be checked by transforming its left-hand side from intrinsic to laboratory variables. This 

10 



yields the relationship for operators: 

u-H^JM = jfiJp - jp D]^{e-\x,)) i^io . (10) 

Upon averaging over eigenfunction "^j of the Hamiltonian H^q ([7]), the second term in the 
right-hand side (1101) yields 0, since \E'/ is separable in functions, dependent on nucleonic and 
additional variables, with the ground-state function of the auxiliary oscillator being real and 
satisfying the equation (/^) = 0. In the end, the expectation value of ffTOj) is equal to the 
right-hand side of ([9]). 

The Hamiltonian ([8]) and operator IJ, involved in the equation (|9]), act in the unified 
space of nucleonic and angular variables, that greatly complicates application of traditional 
methods of the many-body theory. To simplify the problem, the ME of operators /„, acting 
in the space of collective angles (so called geometrical factors), can be found in different ways 
depending on the state spins. The ME for band-head states and close to them rotational 
levels with / ^ i^, can be taken into account exactly. Here, parts of the Hamiltonian, con- 
taining geometrical factors, should be diagonalized in a basis of low dimensionality, which 
is built of states with different values of the quantum number K. The wave functions of 
intrinsic motion in these states are weakly dependent on /, and hence, little differ from wave 
functions of zero-order approximation, i. e. with rotation ignored, so iterative processes in 
self-consistent calculations with account of quasiparticle-rotation interaction will converge 
rather fast (as compared to states with large J). Another, more simple and in many cases 
preferable, way of treating the terms with the operators la directly in the unified space of 
nucleonic and angular variables is to implement a semi-perturbative approximation, consid- 
ered in Sec JIVI For the band-head states and close to them rotational levels, the operator F, 
introduced in this section, is very small that insure the high accuracy of the approximation. 

Good accuracy for rotational states that are not band-head states is achieved with the 
standard SCM approximation; to apply it, one should get rid of collective angles and make 
use of approximate equations: 

4^ V/(/ + i)-(j2), j^^o. (11) 

This is carried out by multiplying the corresponding equations, e. g., those for eigenvalues 
of the Hamiltonian H (or -f^gff), by /^-functions, then integrating the equations over i} and 



summing over K with subsequent use of Eqs. (ITT]) (see Ref. [5^). Note that such approxi 



mation proves in many cases to be acceptable (when calculating energies) even for band-head 

11 



states. Of interest in the present work is mostly the quahtative aspect of the problem, so we 
confine ourselves to the SCM approximation (iTTl) . Then the quantum constraint ([9]) reduces 
to the ordinary SCM condition on the nuclear wave function |$/), dependent on nucleonic 
variables only: 

(J.) = vni + 1) - (J!) ■ (12) 

Consider now the first approximate constraint condition: smallness of the ME of operators 
6a. On the one hand, the smallness arises as a consequence of the limiting process in the 
equation (|H]), resulting in small expectation value {6"^) (and similar averages for higher powers 
of 6). This allows simplifying the initial expressions for the Mikhajlov Hamiltonian and 
commutation relations, and more precisely, their contribution to the variation functional. 
The functional, built on their basis, represents a series in powers of ME 6a and in powers 
of products of operators 6a {h — Jb)- From the conditions of 6 smallness it follows that 
among terms, differing in the powers of ^, only terms of the lowest power contribute to the 
functional. 

On the other hand, to ensure that the condition holds, it is necessary to retain in the 
variation only the terms of the lowest power in 6, i. e. use the expressions (jl]) and ([5]) 
in the functional. Indeed, the ME of operators 6 are determined by minimization of the 
functional {H) on the hypersurface, defined by the equation ([T]). The varying functional is 
a multidimensional polynomial of a very high degree in the ME and can have very many 
minimums. For the condition of 6 smallness to hold, the functional minimum should be 
found, which position is close to zero values of the ME. To find it all polynomial terms of 
power higher than two should be discarded. 

Smallness of 6 does not imply smallness of 6a{Ib — Jb), as: i) by virtue of Eq. f|T3|) . 
coherent sum of 6 a J a over states of all nucleons is of order 1, and ii) the value increases 
with increasing nuclear spin. However, strong dependence of terms of high power in 6 {I — J) 
on spin / is compensated by the same dependence on (J^) with the account of Eq. ( fT2l) . 
Because of the presence of commutators in Eq. (jll), each multiplication by (/ — J) is linked 
with additional coherent summation only in terms of even power in this factor. Here such 
summations appear in expectation values (JaJa)- But the number of couples JaJa is half of 
number of the operators 6, therefore every additional odd power by J gives an additional 
small factor, which can be evaluated as one-particle contribution to {[Ja,i6a]) (see the next 
section). Thus, only three first terms can be retained in the series (II]) for H. 

12 



Finding the operators 9a by the variational method, one, by that very act, define the 
optimal IF (OIF). In our approximation, this is the frame of reference, where a nucleus is 
described to the best advantage by a single-quasiparticle wave function. Similarly, the OIF 
can be defined for a more general form of wave functions (the OIF in the random phase 
approximation for even-even nuclei is briefiy discussed in Ref. 33|); this is, however, beyond 
the scope of the present work. The general scheme of the iterative computation within the 
given approximation is as follows. The ME of operators 9a are obtained by minimizing the 
expectation value {H) under the supplementary condition 

{[Ja,l9a]) = l. (13) 

The initial stage uses eigenfunctions of the SCM Routhian H^^ = H — uJ^ with the angular 
frequency u determined by Eq. f|T2|) . Thus, the variational functional has the form: 

J'={H)+ L{[H,z9a]AJa) + ^-^^^{[[H,^9a],^9a])+Vail-{[Ja,^9t]))\, (14) 

a=x,y ^ ^ 

where rja are Lagrange multipliers, AJ^ = Ja — (Ja)- After determining 9a, the Hamiltonian 
H is constructed, and eigenfunctions of the Hamiltonian if^ = H — uJ^ are sought with the 
condition ( !T2l) . These new functions are used to find more correct 9ai and so on, until the 
convergence of iterative process. 



III. MANY-PARTICLE OPERATORS OF NUCLEAR ORIENTATION ANGLES 

The operators 9a are many-particle ones and so, in the representation of creation and 
annihilation operators a^, a, they have the form 

9 = — 9i2,,,N, l'2'...N' di^l ■ ■ ■ (^N^N' ■ ■ ■ a2'0'V , (15) 

where the ME 9i2,,,n, r2'...N' of the operator 9a are written in the antisymmetrized form, index 
a {a = X, y) is omitted. Direct calculation of the ME by variational method is hardly feasible: 
A^-particle operator in the basis of dimensionality M requires variation of parameters. 
Yet, the situation changes radically, if one takes into account that, with proper selection of 
the basis, only some combinations of the ME actually contribute to the observable quantities. 

A convenient choice of the initial basis in flTSj) is the canonical basis (see jsl, [ssl), obtained 
by the HFB method for Hamiltonian H (hereinafter we consider the Hamiltonian with ef- 
fective forces), or, for the sake of simplicity, single-particle functions, obtained by the HF 
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method for the same Hamiltonian. Then, in representation of quasiparticle operators , a, 
obtained by the HFB method for the SCM Routhian iJ^ (or Hamiltonian H^S)-, operator 6 
has the form 6 = 6^^ + 61^°+°^ + ^40+04+22+31 _|_ _ _ _ _ Hereinafter, for any operator A brought 
to normal form, symbol A^^ denotes its part, containing i operators and j operators 
a. Applications of similar denotation in terms of other creation operators will be always 
accompanied by the pointing to which operators they pertain. 

To derive variational equations for the ME of operators 9a-, it is necessary first to esti- 
mate order of smallness of their different components and use the estimates to simplify the 
functional under variation. We begin consideration with 9]^^'^^ . Comparison of quantities 



{[Ja,i9a]) = 1 and (AJ^) shows that 9^^+^^ ~ Jl^+^^ /2{AJl) on the average, and, since 



less than J^^"*"^*^. In the Hamiltonian (jlj), only each even power of 9 J is linked with coherent 
summation over states of all particles, whereas each odd power is not accompanied by such 
summation due to the presence of commutators and selection rules for the ME of operators 
9 and J. Therefore, terms of power higher than two in 9]^^"^^ J^^^"^^ in (jlj) can be neglected. 

To simplify further, consider the structure of operators 9"^^ In each 9^^ , multipliers before 
normal products of the quasiparticle creation and annihilation operators have the form of 
foldings of the ME 9i2...n,V2'...n' with density matrices pim = {al^ai) and pair tensors Kim = 
(amdi), = (a|aj^); the summands of the operator 9^^ contain products of different numbers 
of matrices p and k. For example, the structure of ^i^+^o is: 



(AJ^) ~ 50 for deformed nuclei, one can conclude that 9^^^"^^ is two orders of magnitude 




Jpp' 



(16) 



X ; dp Cipf '. , 



(17) 



^^ii+20(^t^t + = 1 ^ ij ^^(;v) 9'^^ + ... + gN-2{N) 9. 



■ N-2 



Jpp' 




pp' 



+ 



[g,{N)9'^^+ ... +gN-2{N)9. 



.N-2 




(18) 




(19) 




2...N,2'...N' 



'pp'3...Af,22'3'...Af' ^^22' P3'3 • • • PN'N ■ 



(20) 



14 



Here expressions, like : ap'ttp :, denote normal product of operators, obtained after the Bogoli- 
ubov transformation, fk{N) and gk{N) are combinatorial factors, dependent on N, in par- 
ticular: fi{N) — N, gi{N) — N{N — l)/2. Of similar structure are 6'*-^ with large values 
i+j. 

The less operators and a are contained in component 9^^, i. e. the less i+j, the 
more coherent summations (foldings 0i2...n,i'2'...n' with pim or with Kim) are contained in 
the multiplier of operator . . . a^a ... a. The convohition 9p contains coherent summation 
over all filled nucleonic states. Such sum can be estimated by order of magnitude as one 
term, multiplied by the number of nucleons A. This folding proves to be much greater in 
magnitude than the folding 9k, since the latter contains summation only over single-particle 
states in the vicinity of the Fermi surface where occupation probabilities are smeared out by 
the pairing correlations. Note that with no pairing correlations, one has: 9k = 0. Similar, 
but approximate relation, is useful to simplify the description of highspin states and also 
other nuclear states with slight pairing. 

However, multipliers of operator . . . a^a . . .a contain in addition to the foldings also 
factors fk{N) or gk{N) that describe the number of ways to obtain the foldings. li N = A, 
i. e. the operator 9 is A-particle one and contribution of each nucleon to the operator depends 
on contributions of all other nucleons of the nucleus, then the combinatorial factors outside 
Kim are, in many cases, greater than the multipliers of pi^- Nonetheless, by virtue of the 
condition that 9 are small enough, the contribution of a nucleon to angle of deviation of 
nuclear axes from the intrinsic frame is dependent only on a small number of other nucleons 
and, hence, N < A {as a nice illustration, one can consider an expression for a small angle 
of orientation of the tensor-of- inertia axes for A-particles). For example, in case of the two- 
particle operator 9, all factors fk{N), gk{N) are of the same order of magnitude, and therefore 
terms with ki^ are negligible. There is an analogy here with components of the two-particle 
Hamiltonian: its components, containing Kim, contribute very little to the total energy of the 
nucleus as against those with pi^ (their contributions to excitation energies of the nucleus 
are comparable, but excitation energies of themselves are much less than binding energy). 
Yet, contributions of 9 exhibit an important difference: they are themselves of some order 
of smallness, and the corrections that arise in taking account of components 9, containing 
Kim, a-re corrections to corrections. In a more general case, even if ~ ^4/4, sums with pi^, 
prove to be much greater than those with Kim (see the preceding paragraph) . Such estimate 
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is of qualitative nature. In view of importance of this problem, we will derive a quantitative 
estimate in another independent way. 

To evaluate ^^0+31+22^ differentiate the variational functional with respect to the ME 
of the operators 6^^. Since, by virtue of their structure, ME of 6"^^ are on the average not 
greater than ME of 6*2°+", and the values of ME of 9^^ and 9^^ are on the average of the 
same order of magnitude as the ME of 9^^, it suffices to retain in the functional the terms 
of power not higher than two in these operators. At that (for estimation), there is no need 
to derive all expressions in detail: it suffices to analyze only the structure of the resulting 
linear equation and take into account the validity, on the average, of the following estimate 
for ME: 89^^/39^^ ~ 1 . Then we get that 9^^ ~ 9^^J^^J^y{AJ^). By using the above 
estimate for 6'^°, one can easily make sure that 9^^ << 9"^^ and, hence, operators 9"^^ with 
i+j =4 can be neglected. A more thorough analysis confirms validity of the conclusion even 
in cases, when the average of the transformed Hamiltonian has the number of summations 
over single-particle indices for some terms, containing 6'^°, in excess of that in similar terms, 
containing 9"^^. 

The similar way can be used to derive estimate for 9^^ and show that 9^^ « 9^^. Ex- 
tending the considerations further, we get finally that operators 9"^^ with i + j > 2 can be 
neglected. 

Thus, contribution to observable quantities comes from the operators 

9''+'\a^a) = Y,0ii' (21) 
9''^'\a^a^ + aa) = lj2[ C ■ < ■ + ^ : ] , (22) 

^PP' ^ ^ 9p2...N,p'2'...N' P2'2 ■ ■ ■ PN'N , (23) 

2...N, 2'...N' 

^PP' ^ ^ (^pp'3...N,22'3'...N' f^22' P3'3 • • • pN'N , (24) 

2...N, 2'...N' 

where multipliers outside normal products contain maximum number of density matrices. 

In principle, all coefficients 9^^J,, 9^^, can be found by variational method. This does 
not complicate calculations much as against the approximation, when only operators f l2T]l 
are taken into account. However, the present work will make use of a certain smallness of 
9^'^^ as compared with 9^^^ and neglect the operator fl24|) . To make sure of smallness of 9^'^\ 
compare expressions ( |23l) and (12^ . Although the coefficients 9^^, contain an additional factor 
(A^ — l)/2, the factor is compensated by smallness of the convolution 9k as against the 9p: 
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the former can be approximately estimated as sum over all filled states, and the latter — as 
sum over states near the Fermi surface, participating in pair correlations. But 6^^^, has besides 
an additional smallness, arising for the following reasons. Integration over coordinates in the 
ME 0pp'3,,,N^22'3'...N' binds single-particle wave functions (pp{xi) with (t>2{xi), and functions 
(pp'{x2) with 02'(a;2)- Contribution to expression ( |24|) is made only by states |2), |2') near 
Fermi level, otherwise K22' = 0. But then, for states \p), \p') far from Fermi level, overlaps 
of functions (f)p{xi) with 02(2^1) and (j)pi{x2) with 4'2'{x2) are small and, hence, O'^^, are small. 
Further, the states \p) and |2) must have identical isospin projections and also parity (the 
latter holds because of positive parity of the operator 9). With decreasing pair correlations, 
^22' — ^ and, hence, O'^'^'^ — )■ 0. Under slight mixing of states in the quantum number 
K, defined by the equation Jz\p) = Kp\p), single-particle states in the ME Opp/s...N,22'3'...N' 
should satisfy approximate selection rules: Kp — K2 = ±1, Kpi — K2' = 0. This follows from 
the variational equations for 6, which show that operators 0^'^+'^^ and J^^+^o have identical 
selection rules by K, and this is also demonstrably illustrated by the example of rotation 
angle of tensor-of-inertia principal axes at small values of the angle, when the function 
'arctan' that is involved in the definition of the angle can be replaced by its argument and 
the resulting expression can be rendered in terms of spherical functions. In the model with 
monopole pairing, pair tensor is diagonal: k,22' = ^2'2 ^22 (the bar under index denotes time 
conjugation). Hence it follows that contribution to expression f l24p is made only by two 
states |2), approximately conforming to requirement Kp — K2 = ±1, with the contribution 
suppressed by the multiplier K22 — U2V2, consisting of Bogoliubov transformation coefficients. 
No such restrictions exist for quantities 6'^^^, so, to simplify, one can neglect the operator 
0^^^^*'(a^a^ + aa). It can be mentioned that the K- mixing does not influence essentially to 
the conclusion on the smallness of 6^, because, for a higher mixing, the abovementioned 
arguments may be applied to each component and because the decreasing of the pairing, 

(k) 

caused by the mixing, will diminish 6^^,. 

The quantities 6^^J, are dependent on density matrices and so can change with u, v chang- 
ing. However, within the framework of the semi-perturbative approximation (at not too big 
nuclear spins), used in the next section, the changes can be shown negligible. The operator 
f l2T]) can in such case be considered as some analog of a single-particle operator, i. e. as 
a "quasi-singleparticle" one: '}2pp,{p\0a\p') alap', where {p\6a\p') = {da)ppi and a = x,y. In 
the general case, when the changes are essential, the difficulties, related to such dependence 
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of 6pp,, are easily overcome, if the Hamiltonian H at each iteration is expressed in terms of 
operators a, obtained at the preceding iteration. To clarify this, consider, for example, first 
two steps of the iterative procedure. 

At the first step, after the transformation aj^ = Up^ajj^ + Vp^a^, one obtains the coefficients 
u, V from the HFB equations for the Routhian and then solves the variational equations 
for 9^^J,. At the second step, one uses the same Bogoliubov transformation to express the 
Hamiltonian = H — uJx in terms of aj^, and substitutes the obtained values 6^^,. 
After that, a new transformation ajj = Uf^uO^l + v^yUv of the Hamiltonian H^^ is performed 
with the fixed ME 9^^^, in Eq. (EI]). These M, w -coefficients are determined by HFB equations 
for if^. After the transformation aj, = UpyOl + VpyCty, with the obvious expressions of Upy, 
Vpy via Mp^, Vp^ and u^y, v^y, one obtains the functional in terms of new ME e^P^,, which 
are again determined by the variational equations, but with the new coefficients Upu-, Vpy. 
Thus, such iterative procedure effectively takes into account that the values 9^^^, contain in 
every step not only new ME 6p2,,,N,p'2'...N', but also new density matrices. 

To compare the present method with the approximate AMP (AAMP), based on the 
Kamlah expansion, one can consider some explicit forms for 6 ME and (H). Such expressions 
are given in Ref. [3i|] for ground-state bands of even-even nuclei, but similar formulas may 
be obtained for odd nuclei and excited bands of even ones. For example, consider the 
functional J-" (Eq. (HM ) for the state vector |s) = al|0) with u,v coefficients, satisfying the 
HFB equations for the Routhian (the HFB approximation is not crucial here). Variation 
of J-" with the approximate relations A[Ja,iOh] ~ 0, A[H,i9a] ~ A[Hi^,i9a] yields the ME: 

(0|a, tOa at |0) = (0|a, J, a+ |0) , (25) 

{0\asagapi6aat |0) = ^ (0|a, a, a+ |0) . (26) 

There are the same equations for the Hermitian conjugated ME, but the minus sign appears 
in the right-hand side of the equations. After substituting of these ME, (H) takes the form 
of the well-known AAMP expression for the mean value of the nucleus Hamiltonian (see, 
the review ^). Some more complex consideration makes it possible to obtain, also, the 
AAMP expression without exploiting the Eq. ( ITTj) . But the determination of the quantum 
recoil corrections to wave functions in OIF differs significantly from that of AAMP. Indeed, 
in the next and following iterations of the AAMP method, m, f -coefficients are obtained by 
a minimization of the same expression for the energy (now, with the recoil term), while in 
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the OIF description, u, v are determined by HFB equations for the Hamihonian H^, that 
is, by the minimization of (i/o;). The last functional would coincide with the AAMP one, if 
Eq. ( 125|26|) were valid for all variations of the vector \s). But any variation of \s) changes the 
form of 9 ME. For example, a particular variation + = J2t' ^t' 1^') gives ME {t\i9a\s+6s) 
= Cf {t\i6a\t'), where all ME with t' ^ s are quit different from Eq. fj25| [26]) . To obtain 
the ME with t' s, one has to solve the system of Eq. [2^ with respect to 6*^^, and 
then, with the help of these values, calculate {t\i9a\t'). The difference between AAMP and 
OIF arises, probably, because AAMP is not justified in the cases of strong band mixing. 
More over, it should be mentioned that the variational determination of ^^^1 is much more 
preferable than the determination of the values {Q\ati6a at |0)) (0|«s cxq apiOa o.^ |0): see the 
next section. 



IV. DETAILING OF THE METHOD AND RESULTS OF CALCULATIONS IN A 
SIMPLE MODEL 

As noted in Introduction, main features of the proposed method are conveniently ana- 
lyzed within the framework of a relatively simple, but realistic model, rather than in full-scale 
microscopic calculations. To this end, the present work employs a semi-perturbative approx- 
imation, combining exact diagonalization with perturbation theory. Despite its simplicity, 
such approach has its advantages: the approximation i) allows taking account of the blocking 
effect for states, mixed by the quasiparticle-rotation interaction, ii) ensures orthogonality of 
all resulting states, which poses a certain problem in the presence of terms with Lagrange fac- 
tors Tja and uj in the functional, iii) yields good convergence of iterative procedures and visual 
revealing of roles of various effects, iv) gives the possibility to treat exactly the geometrical 
factors, mentioned in Sec. [TTl Detailed statement of the semi-perturbative approximation 
(with exact account of the geometric factors, but without the recoil operator) along with se- 



lection of parameters of single-particle and pair fields is given in Refs |46|, |56| , so the present 
consideration is confined to description of its features as applied to the OIF method. Most 
equations are given in the framework of the approximations flTTj) since the approximation 
is used in the calculations. The equations for the exact treatment of la operators in the 
semi-perturbative approximation are obtained in close analogy with these ones. 

Consider quasiparticle states with definite values of the angular momentum projection 
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K (with definite values of \K\ and signature), wliicli are obtained by a solution of HFB 
equations for the Hamiltonian if, and operators creating quasiparticles in these states. 
Within the simple model, H contains only a single-particle Hamiltonian, diagonal in the 
representation of the particle operators a|, a^, and pairing interaction, so h\ = uia] + viaj, 
nbl\0) = eib]\0), where n is a quasiparticle Hamiltonian. All the following formulas are 
not limited to the case of this simple Hamiltonian since, with the corresponding replacement 
of Ui,vi by the coefficients of the general Bogoliubov transformation, they are valid for T-L, 
obtained from any energy density functional. 

In the semi-perturbative approximation, eigenfunctions of the Hamiltonian (or, at 
next iterations, H^^) are sought in the form 

^ _giFx/x^Q^t|0) . (27) 
I 

For nucleus states with moderate spin J, the operator F^Ix may be treated as a small 
perturbation. It is useful to chose as one-particle operator, then the transformation 
qiFxIx (ioes not change the expectation value (A^). Expanding the transformed Routhian 
Hi, 

Hi = e-^^^^^ii^e^^^^^ , (28) 

in powers of iF^Ix and neglecting terms of power higher than two, we find operator Fx from 
the equation {H'^)"^^ = that involves only terms with AK = ±1. The terms with AK = 0, 
which do not depend on Fx, become by the definition of amplitudes u,v, while such 
terms, depending on Fx, are small. Hereinafter, the denotations A'-' are used for operators 
A in representation of operators bj. Amplitudes Ci are found by the diagonalization of the 



Hamiltonian {H'^y-^ under the condition ( IT2|) (see equations in Refs. j46|,|56|). As a result, 
the eigenvectors of this Hamiltonian have the form Pj,\0), = ^Ci{p) bj. 

Taking into account these approximations and the smallness of the ME of the operator 
Fx, we get equations of the model. The ME of operator 6*^^ are found by the variation of 
the functional f|T4|) . in which the expectation values are calculated over the state \s) = /3]|0), 
corresponding to the nuclear level in question. Since in the leading term of the functional, 
i. e. in (H), we take into consideration only the single-quasiparticle part of the Hamiltonian 
H, the same approximation is made also in terms of the next order of smallness, i. e., in 
operators, containing 6 and Fx- Then, the variation of the functional ( IT^ leads to the set 
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of linear equations for {l\6J^m) (with I < m) and rj^: 

2 {[2 - - Cl]{ei + em) [LiJ + (Cf - Cl) (e^ - [M-J] tdl 

- J2 { LL (2e™ + + q) + M^^M^Je^ + - 2e„) } zO^ 

- Yl CkCm{ (2q + efc + e„) + M," M,^ (e„ + - 2q) } ^^^f^ 



(29) 



X 

km 



k{kj^l,m) 
1 



((A4)2) 

Here and below the following denotations are used for brevity: 9f^ = (/|6'a;|m), jf^ = {l\jx\m), 
Cl = Ci{s), 



him - 2 I '^JfrnLim ~ Ci^^ CkM^iil^ +Cm ^ CkM^^ili > , 

I A: k ) 

2 "^f,?! — Cl ^ Cfc J^l^ — Cm ^ C'fc 



((AJ.)2) 

- Mm (em + 



Llm. + 



20 
Ik 



(30) 
(31) 



^ ^ Cfc ^jkm + ^ ] C'fc J;^^ — 2Ci Cm {s\Jx\s) 



k k 

The set is appended with equation (JT3l) . which in the said approximation can be written as 

^unmix ^unmix 



0<i<m 



2 E 1^^ 



1 



((AJ.)^) 



20 I 2 rJK 



20 |2 



(32) 
(33) 



0<Z<m 0<Z<m 

Summation in quantities 2'"'^™'^ and jT'^nmix jg Qg^rrjg^j only over those single-particle 
states (with positive signature) that do not take part in the Coriolis mixing of states, occupied 
by an odd nucleon. The "unmixed" are the states that have parity or isospin projection, 
alternative to the corresponding properties of the investigated state \s). For them, Ci{s) = 
Cm{s) = and Eq. f l29l) take a simple form 



i9. 



20 

Im 



1 



(34) 



((AJ,)2) V ei + em, 

which was used in the derivation of Eq. f l32|) . One can included into the set of "unmixed" 
states also non-alternative states, which are not involved in the basis for calculation of the 
amplitudes Ci{s). 
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Similar equations are obtained for and rjy. After solving of all these equation, the 
eigenfunction of the Hamiltonian are sought in the form (1271) . The matrix elements of 
the operator are determined by the equation 

[H, ^Ff ] h - c^tot = , c^tot + {s\[H, i9.M + {s\[[H, iQ^ ^^J k) 4 (35) 

with the obvious solution OF^I^l^ = cototJt^/i^i + Cm)- The amplitudes Ci{s) are found by 
the diagonalization of the Hamiltonian 

(e-^^-^-^,e^^-^-)" = H'' - (utot - {s\[H,i9,]\s) + {s\[[H,i9,],i9,]\s) (slJ.ls))^^ (36) 
+ [H, ^9l'] {s\J^\s) + [H, ^F'^'] 4 - ^ u^.J^' 



2 2v//(/ + l)-(J2) 



The amplitudes are determined self-consistently for each value of Wtot since the Hamiltonian 
itself depends on \s). One can easy check that, after the e^^^^^ transformation, the Eq. ( IT2l) 
depends on Wtot (and on the amplitudes Ci{s)). Therefore it is comfortable, instead of w, 
to find the solution Wtot to this equation. When Wtot and Ci{s) have been obtained, all 
the procedure for 9 and Ci{s) is repeated until the convergence of the iterative process is 
attained. The energy is calculated as {H) (see Eq. [T^ . One can note that, in the given 
model, the Hamiltonian ( l36l) does not contain ME ^f^, while {H) depends on these values. 
Let us make notes about application of the semi-perturbative approximation in the full-scale 
self-consistent calculations with the exact treatment of the operators la- As it was mentioned 
in Sect. [Tll this treatment is usefull in calculations for band-head levels and closed to them 
ones. In this case, the quasiparticle states &[|0) with definite values of \K\ and signature are 
found for the Hamiltonian l-i depending on the operators 9a- The HFB equation for this 
Hamiltonian are obtained by the minimization of the functional ( IT4l) taken without of the 
terms containig the operators la- Then the transformation (1271) is used, where the operator 
Fxlx is changed by ^ Fala and the operators Fa are determined by elimination of the terms 
The Lagrange multiplier ^ in the Hamiltonian = H — laJa is defined by Eq. iQ. 
Return to the simple model. The Hamiltonian is selected as the Nilsson potential with 
constant pairing. In this calculations, we will consider only rotational level energies. At 
present, the theory (particularly, in such simple form) predicts intrinsic motion energies 
with poor accuracy, which can tell on the description of rotational bands both for the worse 
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FIG. 1. Energies of levels of the rotational band 5/2 ^[642] in ^^'^Er. The dashed blue and the dash- 
dot red lines correspond to the results of calculations in ordinary CM and in the semi-perturbative 

n 

version of the OIF method. The solid green line is the experimental data ([57]). 



and for the better. To render the description of rotational states independent of the quality 
of calculations of intrinsic motion energies, these energies are found by fitting of band-head 
energies to experiment, with theoretical values taken for those unavailable from experiment. 
It should be pointed out, however, that all rotational bands of odd nuclei (except for band- 
head states) are calculated without fitted parameters. 

Results of calculations of energies of levels for the rotational band with strong mixing in 
^^^Er are presented in Fig. [H Here, dashed blue, dash-dot red and solid green fines denote, 
respectively, energies obtained without and with account of rotational recoil operator, and 
experimental values ( 57|). As seen from the Figure, taking into account of the effect even 
in the relatively simple version of the proposed model noticeably improves the agreement of 
theory with experiment for the rotational band with strong mixing. The mixing for negative 



parity bands of ^^^Er is less pronounced and so the role of the effect for energies, counted 
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off from l:)aiid-hcad levels, is less than for positive parity ones. For all bands, however, the 
rotational recoil term contributes much to the energy of intrinsic motion; therefore, to take 
correct account of it is of great importance in self-consistent calculations of all quantities in 
deformed nuclei. 



V. CONCLUSIONS 

The present work is an attempt to remove some difficulties arising in the microscopic 
description of deformed nuclei which were mentioned in the Introduction and to develop the 
method of the intrinsic frame which is considered today as too complex and not elaborated 
for calculations. For this purpose, the method of the OIF has been developed on the base 
of the Mikhajlov concept of the unitary transformation to the IF. The method seems to be 
more complex in its justification and development but may be simpler in applications than 
the AMP one, for example, in the self-consistent description of the quasiparticle-rotation 
coupling. Though the obtained results are intended for odd nuclei, most part of them can 
also be useful for even ones. 

Some of deficiencies of the work and ways of feather development are obvious. First of all, 
it is desirable to perform microscopic calculations with effective forces based on contemporary 
energy density functionals. One can hope that such calculations may be applied not only 
for the description of deformed nuclei but also for the study of the effective interactions, 
exploiting for this purpose quasiparticle excitations, and the fact that the single-particle 
motion is much less disturbed by particle in these nuclei than in spherical ones. For the sec- 
ond, the method should be implemented and further developed for calculations of transition 
probabilities. Next, small deviations from axial symmetry can be considered in the frame 
of the given version of OIF by the mixing of states, nevertheless, it is desirable to study 
high asymmetries using the 3D transformation and the corresponding constrains in the IF. 
Also, it should be mentioned that there are many vibrational and quasiparticle-vibrational 
states in the deformed nuclei while the present concept of the OIF is based on the HFB 
method (or on the semi-perturbative approximation) , so the development of the method for 
a self-consistent description of these states is the problem of future. 
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